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ABSTRACT 

We report photometric and spectroscopic observations of the black hole binary GRO 
J1655-40 in complete quiescence. In contrast to the 1995 photometry, the light curves 
from 1996 are almost completely dominated by ellipsoidal modulations from the sec- 
ondary star. Model fits to the light curves, which take into account the temperature 
profile of the accretion disk and eclipse effects, yield an inclination of i — 69.50±0.08 deg 
and a mass ratio of Q — M1/M2 = 2.99 ± 0.08. The precision of our determina- 
tions of i and Q allow us to determine the black hole mass to an accuracy of ~ 4% 
(Mi = 7.02 ± 0.22 Mq). The secondary star's mass is M 2 = 2.34 ± 0.12 Mq. The po- 
sition of the secondary on the Hertzsprung-Russell diagram is consistent with that of a 
« 2.3 Mq star which has evolved off the main sequence and is halfway to the start of the 
giant branch. Using the new spectra we present an improved value of the spectroscopic 
period (P = 2 d 62157±0 d 00015), radial velocity semiamplitude (K = 228.2±2.2 km s^ 1 ), 
and mass function (/(M) = 3.24 ± 0.09 Mq). Based on the new spectra of the source 
and spectra of several MK spectral type standards, we classify the secondary star as 
F3IV to F6IV. Evolutionary models suggest an average mass transfer rate for such a 
system of M2 = 3.4 x 10~ 9 Mq yr _1 = 2.16 x 10 17 g s _1 , which is much larger than the 
average mass transfer rates implied in the other six transient black hole systems, but 
still barely below the critical mass transfer rate required for stability. 

Subject headings: binaries: spectroscopic — black hole physics — X-rays: stars — stars: 
individual (GRO J1655-40) 
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1. Introduction 

A new bright X-ray source (designated GRO J1655- 
40) was discovered July 27, 1994 with the Burst 
and Transient Source Experiment (BATSE) on board 
the Compton Gamma Ray Observatory (Zhang et al. 
1994). The optical counterpart and the radio coun- 
terpart were identified soon after the announcement 
of the discovery by the BATSE team (Bailyn et al. 
1995a; Campbell- Wilson & Hunstead 1994). About 
three weeks after the initial hard X-ray outburst, ra- 
dio jets with apparent superluminal velocities were 
observed emerging from the source (Tingay et al. 
1995; Hjellming & Rupcn 1995), making GRO J1655- 
40 only the second source in our Galaxy observed to 
have superluminal jets and the first one to be optically 
identified. 

Bailyn et al. (1995b) established the spectroscopic 
period (2<?601 ± d 027) and the mass function (3.16 ± 
0.15 Mq) of the system using spectroscopic observa- 
tions made during late April and early May, 1995. 
The photometric observations by Bailyn et al. (1995b) 
from early 1995 showed that GRO J1655-40 is an 
eclipsing binary. The light curves from March and 
April, 1995 all have deep triangular-shaped min- 
ima near the spectroscopic phase 0.75 (caused by an 
eclipse of the disk by the star) and much shallower 
minima near the spectroscopic phase 0.25 (caused by 
an eclipse of the star by the disk). These observations 
confirmed earlier hints that GRO J 1655-40 might be 
an eclipsing binary: Bailyn et al. (1995a) observed 
a single eclipse-like event in the photometry on the 
night of August 17, 1994, and the model of the kine- 
matics of the radio jets proposed by Hjellming & Ru- 
pen (1995) gives an inclination of i = 85dcg. In spite 
of the observations of optical eclipses, there has not 
yet been an unambiguous observation of an X-ray 
eclipse (e.g. Harmon et al. 1995a). The lack of X- 
ray eclipses and the presence of optical eclipses puts 
a tight constraint on the inclination of the orbit. 

Unlike most "X-ray novae," GRO J1655-40 contin- 
ued to have major outburst events in hard X-rays long 
after its initial high-energy outburst (see Harmon et 
al. 1995a). There was an outburst event late March, 
1995 (Wilson et al. 1995) and another one starting 
late July, 1995 (Harmon et al. 1995b). Not surpris- 
ingly, the V magnitude of the source during 1995 was 
typically around V ~ 16.5 (e.g. Bailyn et al. 1995b), 
somewhat higher than the quiescent value of V « 17.3 
(Bailyn et al. 1995a). Furthermore, the shapes of the 



light curves from 1995 are complicated. There are 
night-to-night brightness variations, and the phases of 
some of the optical minima are not aligned precisely 
with the spectroscopic phase (Bailyn et al. 1995b). 
X-ray heating and hot spots on an asymmetric accre- 
tion disk probably play a large role in explaining the 
complex light curves. 

After the July/August, 1995 hard X-ray outburst, 
the source finally settled into true X-ray quiescence. 
From the period of late August, 1995 to late April, 
1996, the source was not detected by BATSE (Robin- 
son et al. 1996). The ASCA X-ray satellite made 
several pointed observations in late March, 1996 as 
part of a large multi-wavelength observing program 
(Robinson et al. 1996). The soft X-ray luminosity was 
found to be quite low, with L x fa 2 x 10 32 ergs s _1 
(assuming a distance of d = 3.2 kpc — see Robinson et 
al. 1996), which is about a factor of 900 times smaller 
than the optical luminosity of the secondary star 
(Z/2 ~ 47 Lq, see Section 6). The extended period 
of X-ray quiescence ended late April, 1996 when the 
all sky monitor on the Rossi X-ray Timing Explorer 
satellite detected a sharp increase in brightness in the 
soft X-rays (Remillard et al. 1996a). The source also 
brightened significantly in the optical and UV (Home 
et al. 1996), and in the radio wavelengths (Hunstead 
& Campbell-Wilson 1996; Hjellming & Rupen 1996). 

We have obtained additional photometry and spec- 
troscopy of GRO J1655-40 February and March, 1996. 
Much of the March data was taken as part of the 
larger multi-wavelength program (Robinson et al. 
1996). These data show that the mean V magnitude 
is consistent with its prc-outburst value and that in 
the optical the system is dominated by light from the 
secondary star. Since the accretion disk contributes 
a small fraction of the light and since GRO J1655-40 
eclipses, we have a unique opportunity to model the 
light curves and obtain a reliable measure of the or- 
bital inclination, thereby leading for the first time a 
reliable mass for a black hole. In Section 2 we describe 
our quiescent photometric observations. We present 
an improved spectroscopic ephemeris in Section 3. In 
Section 4 we discuss the spectral classification of the 
secondary star. In Section 5 we describe models of the 
light curves. Section 6 is the discussion section detail- 
ing the evolutionary status of the secondary star, the 
inclination of the radio jet, and the stability of the 
accretion. A short summary of the paper is presented 
in Section 7. We have also included an Appendix to 
this paper which gives a detailed description of the 
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Table 1 

Journal of Spectroscopic Observations 



1995 April 30-May 4 a ' b 1996 February 24-25 c 



Telescope 
Detector 
Grating 
Resolution 

Wavelength Coverage 
Number of Spectra 



CTIO 4.0 m 
Loral 3072 x 1024 
KPGL #3 (527^/mm) 
3.3 A (FWHM) 
3850-7149 A 
73 



CTIO 1.5 m 
Loral 1200 x 800 
KPGL #3 (527^/mm) 
4.0 A (FWHM) 
4673-6657 A 
12 



a Previously published in Bailyn et al. 1995b. 

b The heliocentric Julian date (HJD) range: 2 449 837.580 to 2 449 841.915 
C HJD range: 2 450 137.792 to 2 450 138.890 



code used to model the light curves. 

2. Observations and Reductions 

The spectroscopic observations from April and 
May, 1995 published in Bailyn et al. (1995b) are sum- 
marized in Table |l|. Additional spectra were taken in 
photometric conditions February 24-25, 1996 with the 
RC spectrograph on the CTIO 1.5 meter telescope. 
The KPGL #3 grating and the Loral 1200 x 800 CCD 
combination gave a dispersion of 1.68 A per pixel. A 
He-Ar lamp was observed repeatedly to give the wave- 
length calibrations. In addition to the 12 spectra of 
GRO J1655-40, we also obtained the spectra of 45 
different stars (most with two or more observations), 
including several radial velocity standards, flux stan- 
dards, and MK spectral type standards (Morgan & 
Kennan 1973). 

Photometry of the source was obtained February, 
March, and April, 1996 with the CTIO 0.9 meter 
telescope and Tek 2048 x 2046 #3 CCD (see Table 
||). All of the nights in February, 1996 were photo- 
metric and about half of the nights in the March- 
April, 1996 run were photometric. Standard IRAF 
tasks were used to process the images to remove the 
electronic bias and to perform flatfielding corrections. 
The programs DAOPHOT He and DAOMASTER 
(Stetson 1987; Stetson, Davis, & Crabtree 1991; Stet- 
son 1992a, b) were used to compute the photometric 
time series of GRO J1655-40 and several field com- 



parison stars. The errors in the instrumental mag- 
nitudes were estimated by computing the standard 
deviations in the light curves of several stable field 
stars. The sizes of these errors were found to be in 
reasonable agreement with the errors reported by the 
DAOMASTER program. We have adopted the errors 
given by DAOMASTER in the analysis presented be- 
low. The instrumental magnitudes were transformed 
to the standard system using previously calibrated 
field stars (Bailyn et al. 1995a). All of the February, 
1996 frames were usable. Out of the all the frames 
taken in March and April, one frame in B, two frames 
in i?, and one frame in / were not used in any of the 
analysis because of cosmic ray hits on or very near 
the image of GRO J1655-40. 

3. A Refined Spectroscopic Ephemeris 

The 73 spectra from April 30-May 4, 1995 were 
rebinned to match the slightly larger dispersion and 
smaller wavelength coverage of the 12 spectra from 
1996. The radial velocities of the 85 spectra were 
found by computing the cross-correlations (Tonry & 
Davis 1979) against a spectrum of 40 Leo (=HD 
89449) , a radial velocity standard star of spectral type 
F6IV (radial velocity = +6.5 ± 0.5 km s _1 ) that was 
observed February, 1996. The cross-correlations were 
computed over the wavelength region between H/3 and 
Ha, excluding the strong interstellar absorption lines 
(see Bailyn et al. 1995a) and regions corrupted by bad 
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Table 2 

Journal of Photometric Observations 



UT Date 


Telescope 


Filter 


Detector 


Number 
of exposures 


1996 February 5,7-20 a 


CTIO 0.9 m 


V 


Tek 2k #3 


24 


1996 February 5,7-20 


CTIO 0.9 m 


I 


Tek 2k #3 


23 


1996 March 21-31 b 


CTIO 0.9 m 


B 


Tek 2k #3 


24 


1996 March 21-31 


CTIO 0.9 m 


V 


Tek 2k #3 


157 


1996 March 21-31 


CTIO 0.9 m 


R 


Tek 2k #3 


24 


1996 March 21-31 


CTIO 0.9 m 


I 


Tek 2k #3 


150 


1996 April l c 


CTIO 0.9 m 


B 


Tek 2k #3 


2 


1996 April 1 


CTIO 0.9 m 


V 


Tek 2k #3 


12 


1996 April 1 


CTIO 0.9 m 


R 


Tek 2k #3 


2 


1996 April 1 


CTIO 0.9 m 


I 


Tek 2k #3 


9 



a HJD range: 2 450 118.836 - 2 450 133.892 
b HJD range: 2 450 163.711 - 2 450 173.895 
C HJD range: 2 450 174.727 - 2 450 174.895 



pixels on the CCDs. Out of the 85 spectra, 84 yielded 
significant cross correlations (V values greater than 
4; see Tonry & Davis 1979). The remaining spectrum 
had a huge cosmic ray hit and was not used in any of 
the analysis presented below. A sinusoid fit, exclud- 
ing the May 2, 1995 data (when the star was partially 
eclipsed by the disk) was performed giving the spec- 
troscopic elements listed in Table |^. The folded radial 
velocities and the best fitting sinusoid are shown in 
Figure [I]. 

We relaxed the assumption of a circular orbit and 
attempted to fit the radial velocities to an eccentric 
orbit. We applied the Lucy & Sweeney test (1971) to 
check the significance of the derived eccentricity (e = 
0.057 ± 0.020). The significance of the eccentric orbit 
fit is much lower than 5%, indicating the eccentric 
orbit fit is no improvement over the circular orbit fit 
to the radial velocities. 

The light curves from February and March, 1996 
folded on the spectroscopic ephemeris are shown in 
Figure |^. The light curves are dominated by ellip- 
soidal variations with maxima at the spectroscopic 
phases 0.0 and 0.5 (the quadrature phases) and min- 



ima of unequal depth at the spectroscopic phases 0.25 
and 0.75 (the conjunction phases). During March, 
1996 we observed three local extrema: a deep min- 
imum on March 22, a maximum on March 24, and 
a shallow minimum on March 26. We fit a parabola 
to the V data from each of these three nights and 
determined the times of the local extrema from the 
fits. The times are (HJD 2,450,160+) 4.835 ± 0.010, 
6.806 ± 0.010, and 8.766 ± 0.010. The spectroscopic 
phases of these three times are 0.261 ± 0.008, 0.013 ± 
0.008, and 0.760 ± 0.008, very close to their expected 
values of 0.25, 0.00, and 0.75. Thus, the phasing of 
the March, 1996 V light curve is consistent with the 
spectroscopic ephemeris given in Table a. 

Our revised spectroscopic period can be compared 
with the period of 2*?2616 ± 0^0016 found by van der 
Hooft et al. (1996) based on photometry from May to 
July, 1995. Our period agrees with theirs to within 
their errors. However van der Hooft et al. (1996) 
then used spectroscopic fiducial points (including the 
one in Bailyn et al. 1995b) and derived the following 
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Table 3 

Orbital Parameters for GRO J1655-40 



parameter result 

Orbital period, spectroscopic (days) 2.62157 ± 0.00015 

K velocity (km s" 1 ) 228.2 ± 2.2 

7 velocity (km s _1 ) -142.4 ± 1.6 

T (spectroscopic) (HJD 2 440 000+) 9 839.0763 ± 0.0055 

Mass function (M©) 3.24 ± 0.09 




-1 -0.5 0.5 

Spectroscopic Phase (P = 2. 63157 days) 



Fig. 1. — The folded radial velocities of GRO J1655-40 and the best fitting sinusoid (see Table ||). Each point has 
been plotted twice. Open diamonds indicate data from April 30, 1995; open triangles, filled stars, and open crosses 
indicate data from May 2, 3, and 4, 1995, respectively. The velocities from February 24, 1996 are indicated by the 
filled triangles, and the filled circles indicate data from February 25, 1996. Parts of this figure appeared in Bailyn 
et al. (1995b). 



ephemeris: 

T min (HJD) = 2449 838.4279(30) + 2.62032(50)7V (1) 

This period differs from our period by 2.5a. Also, the 
phases of the three local extrema in our V light curve 
from March 1996 are by the above ephemeris 0.568 ± 
0.024, 0.320 ± 0.024, and 0.068 ± 0.024, significantly 
different from their expected values of 0.50, 0.25, and 
0.00. Bailyn et al. (1995b) reported great difficulties 
in phasing up their light curves from March and April, 
1995 — no single period could align all of the optical 
minima with the radial velocities. It appears that the 
light curves from May to July, 1995 still suffer from 
phasing problems, although to a much smaller extent 
than the earlier light curves from 1995. The period 



derived solely from the radial velocities appears to be 
the most reliable. 

The refined period of P = 2<?62157 ± d 00015 rep- 
resents a great improvement in accuracy over the pre- 
vious value given in Bailyn et al. (1995b). Using this 
new period, we can now compute accurately the or- 
bital phase of observations made in 1994 and late 
1995. For example, Bailyn et al. (1995a) presented 
photometry from August 17, 1994 which showed evi- 
dence for an eclipse (the source got fainter and redder 
then bluer and brighter over the 4.5 hours of observa- 
tions). Only one such eclipse-like event was observed 
in 1994. Using the refined orbital period, we find 
that the spectroscopic phase of the time of minimum 
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February 1996 



March 1996 




Fig. 2. — The light curves from February, 1996 (left panels) and March, 1996 (right panels) folded on the spectro- 
scopic ephemeris (Table ||) are shown. Each point has been plotted for clarity. Error bars are shown for all of the 
points, but in many cases the size of the errors is smaller than the symbols. 
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Fig. 3. — The rms residuals of a polynomial fit to the difference spectrum as a function of the weight w for three 
different template stars: an F31V star (solid line), an F3V star (dotted line), and a GOV star (dashed line). In each 
case, the minimum of the curve can be accurately computed. The restframe spectrum was from 1996. See the text 
for more details. 



light (estimated to be at the heliocentric Julian date 
2,449,581.63) is 0.80. Curiously enough, the sharp 
optical minimum observed April 2, 1995 (see Bailyn 
et al. 1995b) also has a spectroscopic phase of 0.80 
(the phasing offsets of the light curves will be dis- 
cussed more below). The ASCA X-ray satellite ob- 
served GRO J1655-40 at the end of August 23, 1994 
(Inouc ct al. 1994). During most of this observation 
the X-ray flux stayed relatively low. Also, the X-ray 
spectrum from August 23 was quite different when 
compared to subsequent observations, leading to the 
speculation that an X-ray eclipse was observed (In- 
oue 1995, private communication). It turns out that 
the spectroscopic phases of the August 23 ASCA ob- 
servations are from 0.27 to 0.36 — i.e. the star was be- 
hind the compact object. Thus some other mechanism 
must be invoked to explain the X-ray light curve and 
spectrum from August 23, 1994. 

4. Spectral Classification 

During the April-May, 1995 observing run, the 
spectra of 33 different stars of spectral type F through 
K and luminosity classes V through III were taken. 
During the February, 1996 run, we obtained the 



spectra of 45 different stars, including nine Morgan- 
Keenan (1973) spectral type standard stars of class 
F. We used the technique outlined in Orosz et al. 
(1996) to classify the 1995 and 1996 spectra of GRO 
J 1655-40. First, each continuum-normalized spec- 
trum of GRO J 1655-40 is shifted to zero velocity 
and the lot of them are averaged together, creating 
a "restframe" spectrum. To help minimize the ef- 
fects of an occasional cosmic ray, a "min-max" rejec- 
tion scheme was used where the lowest and highest 
values at each pixel were rejected before the average 
was computed. Then, the velocity shifts between the 
comparison spectra and the restframe spectrum are 
removed. Finally, the steps involved in the compar- 
ison are (1) each continuum- normalized and shifted 
comparison spectrum is scaled by a factor w (where 
< w < 1) and subtracted from the normalized rest- 
frame spectrum, (2) a low order polynomial is fit to 
the difference spectrum, (3) the rms residuals of the 
fit (computed over the same region used to compute 
the cross-correlations) is recorded, and (4) steps (1) 
through (3) are repeated using different scaling fac- 
tors until the minimum rms is found, corresponding 
to the "smoothest" difference spectrum. In finding 
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Fig. 4. — The minimum mis (see the text) is plotted as a function of the spectral type of the template for the four 
different groups. The plotting symbol indicates the luminosity class of the template: open triangles for dwarfs, 
filled circles for subgiants, crosses for giants, and curved squares for supergiants. Templates that had a weight 
to > 1 are not shown (see the text). In each case, the lowest overall rms values occur for the spectral classes F3 
to F5. Note that many of the template stars had two or more observations, so that the total number of template 
stars used is less than the number of symbols shown. 



the rms of the fits to the difference spectra, we were 
careful to avoid interstellar absorption lines and re- 
gions corrupted by bad columns on the CCDs. Fig- 
ure U shows a plot of the rms of the fit to the differ- 
ence spectrum as a function of the weight w for some 
representative cases. In each case, the minimum of 
the curve can be determined accurately. All of the 
above four steps are done on all of the comparison 
spectra, and the comparison spectrum that has the 
lowest overall rms value is taken to be the spectrum 
whose absorption lines most closely match those in 
the restframe spectrum. The relative fluxes of the 
restframe spectrum and the difference spectrum give 
an estimate of the fraction of the total light due to 
the accretion disk. 

We have two restframe spectra of GRO J1655-40, 
one each from 1995 and 1996 (for the 1995 restframe 
spectrum, we did not include the spectra from May 
2, 1995 because the star was partially eclipsed by the 
disk). We also have two sets of template compari- 
son spectra from the same times. Every spectrum in 
each set of templates was compared against both rest- 



frame spectra. Since the 1995 spectra have a slightly 
higher spectral resolution than the 1996 spectra and 
since the 1996 restframe spectrum has a lower signal- 
to-noise ratio than the 1995 restframe spectrum, a 
given template spectrum will have a slightly higher 
minimum rms value when compared against the 1996 
restframe spectrum. We therefore grouped the com- 
parisons into four different cases: the 1995 template 
spectra compared against the 1995 restframe spec- 
trum, the 1995 template spectra compared against 
the 1996 restframe spectrum, etc. The relative values 
of the minimum rms within each group can be used 
to judge the best spectral match to the given rest- 
frame spectrum. Figure |] shows the minimum rms 
value plotted against the spectral type of the template 
spectrum for the four different cases. In all four cases, 
there is a clear trend: the best matches for the GRO 
J1655-40 restframe spectra are the F3-F5 stars. Stars 
with spectral types later than F6 or earlier than F3 
provide much poorer matches to the restframe spec- 
tra. Based on the values of the weight w found for 
the best matches, we conclude that the accretion disk 
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Fig. 5. — The restframe spectrum of GRO J1655-40 from February, 1996 (fourth from the top) and the spectra of 
several F subgiants and giants. The anomalous strengths of the absorption lines near 5900 A and 6300 A in the 
spectrum of GRO J1655-40 are due to interstellar absorption. Each spectrum has been normalized to its continuum 
fit, and offsets have been applied for clarity. 



contributed about 50% of the flux in V during April 
and May, 1995, and less than 10% of the flux in V 
during February, 1996. 

This comparison technique is relatively sensitive to 
the temperature class of the comparison star and rel- 
atively insensitive to the luminosity class of the com- 
parison star. We therefore cannot reliably determine 
the luminosity class of the secondary star from this 
technique alone. However, in the case of GRO J1655- 
40, a main sequence F star is much too small to fill 
the Roche lobe (Bailyn et al. 1995b), so the secondary 
star must be somewhat evolved. Also, in the case of 
the 1996 restframe spectrum, many of the early F 
dwarf template stars were excluded because a weight 
of w > 1 was required to get the minimum value of 
the rms. Hence, we are left with mainly subgiants and 
giants. So, based on these other pieces of evidence, 
we conclude the secondary star in GRO J1655-40 is a 
subgiant, giving its full spectral classification as F3- 
F5 IV. Indeed, with the exception of the strong inter- 
stellar features, it is difficult to tell the difference be- 
tween the 1996 restframe spectrum of GRO J1655-40 
and the spectra of the nearby F subgiants and giants 
(see Figure ||). 



The 1995 spectra of GRO J1655-40 were examined 
in more detail to see if the best fitting spectral type 
depends on the orbital phase. There was no indication 
of any change in the spectral type as a function of 
phase. Further details can be found in Orosz (1996). 

5. Light Curve Models 

The light curves shown in Figure || appear to be 
almost completely ellipsoidal — there are two equal 
maxima per orbit and two minima of unequal depth 
per orbit. The minimum at the spectroscopic phase 
0.25 (when the star is behind the compact object) 
is deeper because the gravity darkening near the L\ 
point is greater and hence the star appears darker (see 
Avni 1978 and the Appendix). The symmetry and 
smoothness of the GRO J1655-40 light curves from 
1996 are in strong contrast to the light curves for the 
other black hole binaries (e.g. McClintock & Remil- 
lard 1986; Wagner et al. 1992; Chevalier & Ilovaisky 
1993; Haswell et al. 1993; Haswell 1996; Remillard et 
al. 1996b; Orosz et al. 1996), where the light curves 
are complicated by "flickering" about the mean light 
curve and by large asymmetries in the light curve that 
slowly change with time (e.g. Haswell 1996). 
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Table 4 
Input Model Parameters 



parameter description comment 



i 


inclination 


free 




Q 


mass ratio (M1/M2 > 1) 


free 






outer radius of the disk a 


free 




'inner 


inner radius of the disk a 


fixed at 


0.005 


?disk 


disk temperature at the outer edge 


free 




t 


power law exponent on disk 
temperature radial distribution 


free 




Prim 


flaring angle of the disk rim 


free 




L.i; 


X-ray luminosity of the compact object 


fixed at 


0.0 


w 


X-ray albedo 


fixed at 


0.5 


T e ff 


polar temperature of the secondary 


fixed at 


6500 K 




gravity darkening exponent 


fixed at 


0.25 


u(X) 


linearized limb darkening parameter 


fixed b 




P 


orbital period 


fixed at 


2 d 62157 


f(M) 


mass function 


fixed at 


3.24 M 


n 


ratio of rotational angular velocity 






to orbital angular velocity 


fixed at 


1.0 


f 


Roche lobe filling factor 


fixed at 


1.0 



a The disk radii are scaled to the effective Roche lobe radius of the compact 
object. 

b Limb darkening parameters from Al-Namity (1978) and Wade & Rucin- 
ski (1985). 



Because the light curves from February and March, 
1996 are smooth and symmetric and because the lu- 
minosity of the star dominates (i.e. the fraction of the 
light from the disk is small, see Section we have 
a unique opportunity to model the light curves and 
obtain a reliable constraint on the orbital inclination. 
We have developed a detailed code based on the work 
of Avni (Avni & Bahcall 1975; Avni 1978) to model 
the light curves, which is fully described in the Ap- 
pendix. The code uses full Roche geometry to account 
for the distorted secondary, and light from a circular 
accretion disk is included. The code also handles mu- 
tual eclipses by the star and the disk. The effects 
of X-ray heating on the secondary star are included 
as well. There are several input model parameters 
which are summarized in Table ^. See the Appendix 
for detailed discussions of these parameters and their 



meaning. 

For simplicity, we have fixed several of the input 
parameters at reasonable values. For example, we 
know the spectral type of the star (see the previous 
section), from which we can get its effective temper- 
ature. The effective temperature of 6500 K is appro- 
priate for an F5 IV star (Straizys & Kuriliene 1981) 
and we will adopt the effective temperature of 6500 
K as the polar temperature. The gravity darkening 
exponent is fixed at 0.25 since the star has a radia- 
tive envelope (see the Appendix) . The limb darkening 
coefficient u(A) is better determined from model at- 
mosphere computations, and we used values interpo- 
lated from tables given by Wade & Rucinski (1985). 
For fits to the 1996 light curves, we have fixed the 
X-ray luminosity at 0, based on the ASCA observa- 
tion in early 1996 that found L x ss 2 x 10 32 ergs s _1 , 
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Table 5 

Fits to the March, 1996 photometry 



parameter 


model 1 


model 2 a 


model 3 


i (degrees) 


69.50 ±0.08 


74.7 ± 1.2 


69.54 ±0.08 


Q 


2.99 ±0.08 


3.3 ±0.5 


3.50 ±0.08 




0.747 ±0.010 


0.71 ±0.01 


0.748 ±0.011 


Tdisk (K) 


4317 ± 75 


4525 ± 102 


4429 ± 75 


Am (degrees) 


2.23 ±0.18 


1.54 ±0.4 


2.91 ±0.18 


£ 


-0.12 ±0.01 


-0.22 ±0.01 


-0.13 ±0.01 




0.0 (fixed) 


0.0 (fixed) 


0.0 (fixed) 


w 


0.5 (fixed) 


0.5 (fixed) 


0.5 (fixed) 


Toff (K) 


6500 (fixed) 


6500 (fixed) 


7000 (fixed) 


X. u 


1.1551 


3.2391 


1.1816 



No checking for eclipses done. 



almost three orders of magnitude lower than the op- 
tical luminosity of the secondary star {Li « 47 Lq, 
see the next Section). In the case where L x = 0, the 
exact value of the X-ray albedo W is of course irrel- 
evant. For fits to the March 1995 light curves where 
L x 3> L2, we have adopted W = 0.5 (see the Ap- 
pendix), although other values of the X-ray albedo 
W have been used (for example van der Hooft et al. 
(1996) used W = 0.4). Since there is mass transfer 
taking place, we assume the star fully fills its Roche 
lobe and that it is in synchronous rotation. 

For fits to the 1996 light curves, we are left with 
six free parameters. The six free parameters in the 
model are the inclination of the orbit (i), the mass 
ratio (Q = M1/M2), the outer radius of the disk 
in terms of the primary's effective Roche lobe ra- 
dius (r<j), the temperature at the outer edge of the 
disk (Tdisk), the flaring angle of the rim of the disk 
(Aim, where the thickness of the disk at the outer 
edge is given by 2^ tan Aim), and the power law ex- 
ponent on the disk's radial temperature distribution 
(£, where the disk temperature varies with radius as 
T{r)=T d ^{r/r d )t). 

We first fit models to the BVRI light curves from 
March, 1996, as the phase coverage and spectral cov- 
erage from that time is the most complete. We folded 
the data on the photometric phase convention used in 
the code where the photometric phase 0.0 corresponds 



to the time the secondary is directly in front of the 
compact object (i.e. Xb(photo) = T (spect) +0.75P). 
As discussed in Section 3, the phases of the three 
local extrema observed in March, 1996 are all late 
in phase by w 0.011. We found that the model fits 
were slightly better after a phase shift of 0.011 was 
removed from the folded curves. Since there are far 
fewer points in the B and R filters, we gave each point 
in B and R seven times more weight when comput- 
ing the chi-square of the fit. This procedure resulted 
in equal weight being given to each filter. Figure || 
shows the fits and the residuals and Table || lists the 
model parameters (under "model 1"). Note that the 
data in all four filters were fit simultaneously. The 
model fits the observed light curves from all four fil- 
ters quite well — the scatter in the residuals is typically 
less than 0.02 magnitudes. The la statistical errors 
for each parameter were estimated using primarily a 
Monte Carlo "bootstrap" method (Press et al. 1992 — 
see the Appendix). For this model, the accretion disk 
contributes about 5% of the flux at 5500 A, consistent 
with the measurement of < 10% from late February, 
1996 (see the previous Section). 

As a test of the accuracy of the la errors derived 
from the bootstrap analysis, we performed fits the the 
March 1996 light curves with the value of Q at fixed 
at several different input values. The other five pa- 
rameters were free and adjusted to find the minimum 
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. 6. — The model fits to the B, V, R, and I light curves from March, 1996 (large panels) and the residuals of 
fits (i.e. the data minus the model — small panels). Table || gives the model parameters under "model 1." 
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Fig. 7. — Upper left: The value of the reduced chi-square as a function of the mass ratio Q, where the value of Q is 
fixed at each input value and the other five parameters are adjusted until the chi-square is minimized. Lower left: 
An expanded view of the xl vs - Q curve. The dotted lines indicate the change in the chi-square required for a 1, 
2, and 3a change in a single parameter (i.e. A% 2 = 1, 4, and 9). Upper and lower right: similar to the left panels, 
but for fits with the inclination angle i fixed at various input values. 



reduced chi-square at each input value of Q. The left 
panels of Figure shows the xl vs - Q curve. The 
dashed lines shown in the lower left panel of Figure 
^ indicate the change in xl representing a 1, 2, and 
3cr change in a single parameter (i.e. Ax = 1, 3, and 
9). A change in Q of 0.08 (the la error derived from 
a Monte Carlo bootstrap analysis) in either direction 
from its value at the minimum causes the xl to in- 
crease by approximately the amount indicated by the 
lowest dashed line. The right panels of Figure ^ show 
the equivalent computation done with the value of i 
fixed at several different values. Like before, a la 
change in the parameter i (0.08) from the best value 
causes the value of xl to increase by approximately 
the amount indicated the lowest dashed line in the 
lower right panel. The xl v - s - * curve has a secondary 
minimum at i = 67.7 deg and xl = 1-1620, which is 
greater than the 2a deviation from the primary min- 
imum at i = 69.5 deg and xl = 1.1551. While we 
are confident that the errors derived from the boot- 
strap analysis are a reasonable representation of the 
true internal statistical errors, there may be system- 
atic errors due to physical effects not included in the 
model. Given how well the model fits the data, we 



suspect that such systematic errors are likely to be 
small. We are also confident that we have found the 
global minimum of xl since we have searched a large 
amount of parameter space. 

For each model with Q fixed at a given input value, 
there exists a value of the inclination i found from the 
best fit. For the whole range of Q shown in the upper 
left panel of Figure ^, the corresponding best values 
of i range from 68.40 < i < 71.1 deg. Likewise, for 
the models with the inclination angle i fixed at sev- 
eral values (the upper right panel of Figure [?]), the 
corresponding best values of the mass ratio Q are in 
the range 2.52 < Q < 4.23. Thus, the key geomet- 
rical parameters i and Q vary little and seem to be 
reasonably well determined. Based on the behavior of 
X 2 near the minimum, we adopt the la and 3a ranges 
of the inclination i: 



69.42 deg < i 
69.0 deg < % 

and the mass ratio Q: 

2.91 < Q 
2.60 < Q 



< 69. 58 deg (la) 

< 70.6 deg (3a) 

< 3.07 (la) 

< 3.45 (3a). 
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Fig. 8. — The system as it appears projected onto the plane of the sky for the model parameters given in Table 
model 1. The photometric phase is Odeg for the bottom panel and 170 deg for the top panel. 
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The geometry of the system is shown schematically 
for two different phases in Figure |[ Note that the 
compact object is not eclipsed by the secondary star. 
We would therefore not expect to see X-ray eclipses 
(assuming the X-rays are emitted from a relatively 
small region centered on the compact object), and in 
fact no X-ray eclipses have been seen (e.g. Harmon et 
al. 1995a). The star does block part of the disk (and 
vice- versa) at certain phases, so we expect to still see 
grazing optical eclipses, the depths of which depend 
on the size of the disk and the relative brightnesses 
of each component. In order to assess how important 
the grazing eclipses are in determining the final overall 
shapes of the light curves, we performed a model fit to 
the BVRI data from March, 1996 where the eclipse 
checking routines were turned off. In this case the 
model light curve is just the normal ellipsoidal light 
curve from a Roche lobe filling star plus a constant 
flux from the accretion disk (where the amount of disk 
light relative to the amount of star light in general 
depends on the filter bandpass). Table || gives the 
values of the parameters (under "model 2"). The fit is 
poor, as one can deduce from the relatively large value 
of the reduced chi-square (xf, = 3.239, compared to 
xl = 1.1551 for model 1). The fits are poor when 
eclipses are not accounted for in the computations for 
two reasons: (i) the difference in the depths of the 
two minima in the model light curve are too small, 
and (ii) the dependence of the light curve amplitude 
with color is not properly accounted for by the model 
(the amplitude of the observed light curve is largest 
in B and smallest in /, see Figure ||) . 

Our numerical experiments also show that the fi- 
nal fitted parameters do not depend strongly on the 
assumed secondary star parameters. For example, in 
the case of a Roche lobe-filling star, the temperature 
will not be constant over the surface (see Equation 
(A3) of the Appendix). In particular, the side of the 
star near the L\ point should be the coolest. The 
temperature of 6500 K corresponding to the observed 
spectral type therefore represents some kind of flux- 
weighted average of the temperatures over the sur- 
face, and technically not the star's polar tempera- 
ture. However, our spectral coverage was not uni- 
form in phase and spectra used to determine the spec- 
tral type mainly came from phases near the quadra- 
ture phases. The temperatures corresponding to the 
brightest parts of the projected stellar disk at these 
phases are within w 200 K of the polar temperature. 
Also, we could not detect any change in the spectral 



type of the individual observations as a function of 
phase (this includes the data from May 2, 1995 when 
the star was behind the disk and the L\ point was 
in full view — see Orosz (1996)), which suggests other 
factors are at work when the observed spectrum at 
a given phase is produced. Thus, we believe we are 
not making a large error when adopting the temper- 
ature derived from the observed spectral type as the 
polar temperature. In this same vein, we fit a model 
to the March, 1996 data where the polar tempera- 
ture of the secondary star was 7000 K, appropriate 
for an F2 IV star (Straizys & Kuriliene 1981). The 
results for the geometrical parameters, given in Table 
H under "model 3" , are quite similar to the results 
using T ff = 6500 K (model 1). Finally, we compared 
the linearized limb darkening parameters u(X) inter- 
polated from tables given by Al-Naimiy (1978) and 
Wade & Rucinski (1985). The shapes of the u(\) 
curves were very similar — the only large differences 
occurred for wavelengths blueward of the B band. As 
a result, the fits to the BVRI light curves of GRO 
J 1655-40 were similar when the limb darkening coef- 
ficients were used from either source. 

The stellar evolution models (see the Discussion) 
indicate the secondary has a radiative envelope. We 
have therefore adopted a gravity darkening exponent 
of (3 — 0.25. We tried a model fit to the light curves 
with the gravity darkening exponent set to (3 — 0.08, 
the value appropriate for stars with convective en- 
velopes (Lucy 1967). The model fits were much worse 
in this case (we could not find a fit with xt < 2). 
When (3 — 0.08, the temperature contrast over the 



star is smaller (see Equation (A3) of the Appendix), 
leading directly to a smaller brightness contrast, es- 
pecially between the side of the star near the L\ point 
and the opposite side. Since the brightness contrast 
between the side facing the compact object and the 
opposite side is greatly reduced, the inclination needs 
to be higher to produce the same observed ampli- 
tude of the light curves. However, as the inclination 
grows closer to 90 deg, the eclipses become deeper and 
the fine balance between the primary and secondary 
eclipse depths is disturbed, making the simultaneous 
fits to all four colors much more difficult. Thus, based 
both on the stellar evolution models and the quality 
of the fits, we conclude that the gravity darkening 
exponent of (3 = 0.25, appropriate for a star with a 
radiative envelope, is the correct choice. 

As one can see from Figure 0, the shapes of the 
February, 1996 light curves are almost identical to 
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Fig. 9. — The two panels show the February, 1996 data minus the March, 1996 model (Table [|, model 1). 



the shapes of the light curves from March, 1996. The 
only noticeable difference between the two curves oc- 
curs near the minimum at the photometric phase 0.5 
(which is the spectroscopic phase 0.25): that mini- 
mum in the February light curves is slightly deeper 
than that minimum in the March light curves. Not 
surprisingly, models that fit the March, 1996 data are 
consistent with the February, 1996 data (Figure |^): 
the residuals in the V band have a mean of 0.009 and 
a standard deviation of 0.019 and the residuals in the 
I band have a mean of 0.013 and a standard deviation 
of 0.013. The only deviant points in the two filters are 
those near the photometric phase 0.5 which are about 
0.05 magnitudes fainter than the model that fits the 
data from March. 

The light curves of GRO J1655-40 from March and 
April, 1995 (see Bailyn et al. 1995b) look quite differ- 
ent from the light curves from 1996 (Figure |J). First 
of all, the source was about 0.7 to 0.8 magnitudes 
brighter in V during early 1995 than it was during 
early 1996, probably a direct result of the fact that 
there was still considerable hard X-ray activity tak- 
ing place during that time (e.g. Wilson et al. 1995). 
In 1995, the minimum near the spectroscopic phase 
0.75 (when the secondary star is in front of the com- 



pact object) is the deeper of the two, opposite of the 
situation in the 1996 light curves. There are large 
asymmetries in the 1995 light curves, and the phases 
of the minima of the three light curves shown in Bai- 
lyn et al. (1995b) do not all line up properly with 
the spectroscopic phase (the difference in phase is 
as large as 0.05). Similar (but smaller) phase shifts 
in the optical light curve minima have been seen in 
some low mass X-ray binaries (for example X1822- 
371, Hellier & Mason 1989). The phase shifts of the 
minima are presumably caused by large distortions in 
the accretion disk. Hellier & Mason (1989) invoked 
a thick accretion disk where the thickness of the rim 
depended strongly on the azimuth angle to explain 
the X-ray and optical light curves of X1822-371. It 
is quite clear that our simple model where the disk 
is a flattened azimuthally symmetric cylinder will not 
be able to produce a model light curve whose minima 
are shifted from their expected phases. 

Out of the three light curves from 1995 presented 
in Bailyn et al. (1995b), the V and I light curves 
from March 18-25, 1995 are the most symmetric. In 
addition, the minima of these light curves line up ap- 
proximately at their expected phases when the data 
arc folded on the revised spectroscopic ephemeris prc- 
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Photometric Phase 

Fig. 10. — The models fits to the V and / light curves from March, 1995 (top panels) and the residuals of the fits 
(i.e. the data minus the model — bottom panels). Table ^ gives the model parameters. Note the scale changes in 
the plots compared to Figure @. These data were published previously in Bailyn et al. 1995b. 



sented in Table || (whereas the minima of the later 
light curves do not). We attempted, a model fit to the 
March 18-25, 1995 data to see if the overall shape of 
those two light curves can be explained by our model 
(as discussed above we will not be able to model all of 
the 1995 light curves because of the large phase shifts 
of the minima of the other light curves). For fits to 
the March, 1995 light curves, we must take into ac- 
count X-ray heating of the secondary star. The part 
of the secondary star facing the compact object will 
be hotter and hence brighter than it would otherwise 
be in the absence of any external heating. We set the 
X-ray albedo at W = 0.5 and let the X-ray luminos- 
ity of the compact object L x be a free parameter. We 
fixed the inclination at i = 69.50 deg and the mass 
ratio at Q — 2.99, the values found from the model 
fit to the March, 1996 light curves. The results are 
listed in Table ^, and Figure [lC] shows the fits. 

Considering the simplicity of the model, the over- 
all shapes of the V and I light curves are fit reason- 
ably well. There is considerable "flickering" about the 
mean light curve with deviations larger than the ob- 
servational errors, resulting in the large value of the 



reduced chi-square. Note the difference in the disk 
parameters between the fits to the March, 1995 and 
March, 1996 data: the disk is much larger and hotter 
in 1995. It is interesting to note that the disk and 
the star contribute roughly equal amounts of flux in 
the V band for the model that fits the March, 1995 
data, which is in good agreement with the value of the 
disk fraction estimated from the spectra from April 
and May, 1995 (which were not simultaneous with the 
photometry). The observed X-ray luminosity, as de- 
termined from the BATSE daily averages, varied by 
about a factor of four between March 18 and March 



25, 1995, with 5.7 x 10 36 £ L x £ 2.3 x 10 37 



erg s 



(S. N. Zhang, private communication). Considering 
the simplistic way the X-ray heating is handled in the 
model, the fitted value of L x — 3.7 x 10 36 ergs s _1 is 
remarkably close to the observed range. A slightly 
lower value of the X-ray albedo W (like the value of 
W = 0.4 used by van der Hooft et al. 1996) would re- 
quire a larger value of L x t o pro duce the same heating 
effect (e.g. see Equation ( A23 ) of the Appendix), so 
we do not consider the slight difference between the 
observed X-ray luminosity and the fitted luminosity 
to be a problem. Finally, we note that the depth and 
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Table 6 

Fits to the March, 1995 photometry 



parameter 


March, 1995 
model 4 


i (degrees) 


69.50 (fixed) 


Q 


2.99 (fixed) 




0.95 ±0.01 


Tdisk (K) 


6694 ± 100 


llm (degrees) 


12.6 ±0.1 




-0.18 ±0.01 


L x (ergs s x ) 


3.72 x 10 36 


W 


0.5 (fixed) 


T cB (K) 


6500 (fixed) 


xl 


66.0 



shape of the minimum at the photometric phase 0.5 is 
sensitive to the amount of X-ray heating. Given that 
the observed hard X-ray flux from the source varied 
by about a factor of four during March 18-25, 1995, 
we should not be surprised to see a large amount of 
scatter in the optical light curves around phase 0.5. 
In summary, although we cannot fit every detail in 
the March, 1995 light curves, we conclude that the 
basic overall shapes of the light curves (i.e. the rela- 
tive depths of the minima) can be explained by the 
effects of X-ray heating on the secondary star. To ex- 
plain the phase-shifted minima of the other 1995 light 
curves, we would probably have to invoke an accretion 
disk whose shape and brightness profile are grossly 
distorted — something which is beyond the scope of 
our current model. 

The code also computes radial velocities. For a 
Roche-lobe filling star, the center of light does not 
precisely follow the center of mass of the star, result- 
ing in a slightly distorted radial velocity curve (for 
example, see Kopal 1959 or Wilson & Sofia 1976). 
Partial eclipses of the star can result in an asym- 
metric contribution to the star's rotational velocity 
to the overall observed radial velocity, leading to a 
distorted radial velocity curve near the conjunction 
phases. On May 2, 1995, the secondary star was ob- 
served spectroscopically by Bailyn et al. (1995b) as it 
moved through the spectroscopic phase 0.25 (i.e. as 
it moved directly behind the disk). The radial veloci- 
ties from that night (Figure [l]) deviate systematically 



from the sine curve: just before the phase 0.25, the 
radial velocities are larger than expected and just af- 
ter the phase 0.25 the radial velocities are smaller 
than expected. Stated another way, the residuals of 
the sine curve fit to the radial velocities (in the sense 
of "data-model") are positive just before phase 0.25 
and negative just after phase 0.25. As we discussed 
above, the accretion disk must have had a complex 
and variable shape to account for some of the features 
of the 1995 light curves. Modelling the radial velocity 
curve of GRO J1655-40 is made more difficult by the 
lack of simultaneous photometry. Nevertheless, the 
model shown in Figure [l^ does produce radial veloc- 
ity residuals with the correct shape near phase 0.25 
(positive just before phase 0.25 and negative just af- 
ter), but the amplitude of the model residuals was 
only about 15 km s _1 , whereas the observed resid- 
uals have a much larger amplitude of « 90 km s _1 . 
Indeed, we could not explain the large amplitude of 
the radial velocity residuals from May 2, 1995. Fur- 
ther details of the radial velocity curve models can be 
found in Orosz (1996). 

6. Discussion 

The precision of our results for i and Q allow us 
to determine the masses of the components to un- 
precedented accurately for a black hole system (see 
Table fj]). Indeed the formal precision of our results 
is so high that we expect the true uncertainty of our 
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Table 7 

Component masses for GRO J1655-40 



confidence 
limit 


i 

(degrees) 


Q 


Mi 
(M Q ) 


M 2 
(Mq) 


la 

3(7 


69.42-69.58 
69.00-70.60 


2.91-3.07 
2.60-3.45 


6.80-7.24 
6.42-7.63 


2.20-2.46 
1.86-2.90 



results comes from physical effects not included in 
the models. However, given the good fit between our 
models and the data, we expect that the systematic 
errors are also quite small. In the following we discuss 
the implications of our results for the nature of the 
secondary star, the radio jet, the evolutionary history 
of the system, and the outburst mechanism. We will 
use the parameters listed in Table ^ for the system 
throughout the following discussion. 

6.1. Luminosity and Nature of the Secondary 
Star 

The mass ratio of Q = 2.99 ± 0.08 we found for 
GRO J 1655-40 is the smallest among the seven black 
hole systems with low mass secondaries. Much larger 
values of Q have been found for the three systems 
where the mass ratio has been directly computed via 
the measurement of the rotational broadening of the 
secondary star's absorption lines: Q — 14.9 for both 
A0620-00 (Marsh, Robinson, & Wood 1994) and V404 
Cyg (Casares 1995), and Q = 23.8 for GS 2000+25 
(Harlaftis, Horne, & Filippenko 1996). There are 
hints that the mass ratio may be extreme in XN 
Mus91 (Orosz et al. 1994), and GRO J0422+32 (Fil- 
ippenko, Matheson, & Ho 1995) as well. In GRO 
J1655-40, we have a way to independently check the 
value of the mass ratio by comparing the observed 
luminosity of the secondary star to the luminosity in- 
ferred from the model fits. To compute the observed 
luminosity in V, we adopt a mean V magnitude of 
V — 17.12, which is the mean of the model fit to the 
March, 1996 data, a distance of d = 3.2 ± 0.2 kpc 
(which is tightly constrained by the kinematics of the 
radio jets — see Hjellming & Rupen 1995), a disk frac- 
tion in V of 5% ±2% (Section 4), and a color excess of 
E(B-V) = 1.3±0.1 (Horne et al. 1996; Horne private 
communication 1996). Assuming a visual extinction 
of A v = 3.1E(B — V), the observed luminosity in V of 



the GRO J1655-40 secondary is L ohs = 46.6±13.6 L Q . 

The luminosity of the secondary star inferred from 
the models is easily derived from some simple rela- 
tions. Given the measured masses and orbital period 
of the system, we can use Kepler's Third Law to deter- 
mine the semi-major axis. Eggleton's (1983) expres- 
sion for the effective radius of the Roche lobe then 
determines the size of the secondary, and the temper- 
ature inferred from the spectral type and the Stefan- 
Boltzmann Law determines the luminosity. Figure 11 
shows the computed mass, radius, and luminosity of 
the secondary star as a function of Q, assuming an in- 
clination of i — 69.5 deg, an effective temperature of 
T eS = 6500 K, and an orbital period of 2?62157. The 
vertical dashed line in the three panels indicates the 
value of Q from the best model fit, and the horizontal 
dashed and dotted lines in the lower panel indicate the 
observed luminosity of the GRO J1655-40 secondary 
and its ler error. Evidently, models with Q <^ 3.5 are 
needed to produce a secondary star with a luminosity 
large enough to match the observed luminosity of the 
GRO J 1655-40 secondary. Models with larger mass 
ratios (i.e. Q J> 5) not only produce under- luminous 
secondary stars, but also provide fits to the data with 
relatively large xl values (Figure [?]). The agreement 
between the luminosity of the secondary calculated 
from the orbital parameters and that inferred from 
the known distance and reddening of the source pro- 
vides a satisfying consistency check for our models. 

The observed luminosity of the GRO J1655-40 sec- 
ondary star is a fairly well-constrained quantity, ow- 
ing to the fact that the distance is well determined 
from the kinematics of the radio jets (Hjellming & 
Rupen 1995) and the amount of interstellar extinc- 
tion is well determined from high quality UV spec- 
tra obtained with the Hubble Space Telescope (Horne 
et al. 1996). The observed spectral type of the star 
provides a fairly good measure of its effective tem- 
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Fig. 11. — The mass (top), radius (middle), and the luminosity (bottom) ol the secondary star as a lunction of the 
mass ratio Q, assuming T cS = 6500 K, i = 69.5 deg, f{M) = 3.24 Mq, and P = 2 d 62157. The dashed and dotted 
lines in the lower panel indicate the observed luminosity of the secondary star in GRO J1655-40 and its la error. 



perature (e.g. Straizys & Kuriliene 1981), and we 
adopt T e g = 6500 ± 250 K. We therefore can accu- 
rately place the secondary star of GRO J1655-40 on a 
Hertzsprung- Russell diagram (Figure |l2|). Using the 
Yale Stellar Evolution Code (Guenther et a!. 1992) 
with updated opacities (Iglesias & Rogers 1996), we 
computed the evolutionary tracks of a 2.1 Mq star, 
a 2.3 Mq star, and a 2.5 Mq star, assuming a so- 
lar metallicity, a hydrogen abundance of 71%, and a 
mixing length coefficient of 1.7, all typical values for 
population I stars. These evolutionary tracks are also 
shown in Figure [l2|. The secondary star falls very near 
the track of the 2.3 Mq star, roughly halfway between 
the main sequence and the giant branch. This pro- 
vides yet another consistency check for our model. It 
is interesting to note the radius of the 2.3 Mq star 
near the position of the secondary: 690 Myr after the 
zero-age main sequence (ZAMS) the star has a ra- 
dius of 4.85 Rq (and a temperature of T e g = 6813 K 
and a luminosity of L 2 = 45.8 Lq); 692 Myr after 
the ZAMS, the star has a radius of 5.21 Rq (and 
T cff = 6500 K and L 2 = 43.8 Lq). For Q = 2.99±0.08 
and i = 69.50±0.08 deg, the effective radius of the sec- 
ondary star is 4.85 ± 0.08 Rq (when given the values 



of the orbital period and mass function listed in Table 
^). Thus, the observed luminosity, the observed tem- 
perature, and the inferred radius of the GRO J1655- 
40 secondary star are all consistent with a ss 2.3 Mq 
single star that has evolved « 690 x 10 6 years past 
the ZAMS. 

6.2. Inclination of the Radio Jet 

Hjellming & Rupen (1995) showed that the rel- 
ativistic radio jets GRO J1655-40 possessed during 
late 1994 were inclined 85 ± 2 deg to the line of sight. 
In addition, the jets were apparently rotating about 
the jet axis at an angle of 2 deg and with a period 
of 3^0 ± 0^2. Our light curve models show that the 
orbital plane is inclined 69.5 deg to the line of sight. 
Hence, the jet axis was tilted by 15.5 deg from the 
normal of the orbital plane. If the jet axis was per- 
pendicular to the orbital plane, the 2 deg offset might 
be explained by a warped accretion disk near the com- 
pact object. In this case, the roughly three day peri- 
odicity would be closely related to the orbital period 
of the binary since the orientation of the warp would 
depend on the orbital phase. However, since the jet 
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Fig. 12. — A Hertzsprung-Russcll diagram showing the position of the secondary star of GRO J1655-40. Evolu- 
tionary tracks are shown for a 2.1 Mq star (lower curve), a 2.3 Mq star (middle curve), and a 2.5 Mq star (upper 
curve), assuming a solar metallicity. In each case, core hydrogen exhaustion takes place just after the tiny loop. 
The computations of the tracks for the 2.1 Mq star and 2.5 Mq star were stopped short of the giant branch. 



axis is not perpendicular to the orbital plane, it ap- 
pears difficult to account for both the 15.5 deg tilt in 
the jet axis and the 2 deg rotation about the axis by 
invoking a simple warp in the inner accretion disk. If 
the inner disk were warped by 15.5 deg and the axis of 
the jet were perpendicular to the warp, we might ex- 
pect to observe a large precession of jets because the 
orientation of the warp depends on the orbital phase. 

Alternatively, the 15.5 deg tilt of the jet axis may 
be a measure of the inclination of the spin axis of 
the black hole with respect to the orbital plane. In 
this hypothesis the jet would need to be collimated 
in some way by the Kerr geometry of the black hole 
itself. The "wiggles" observed in the jets might be due 
to the low mass ratio of the system which results in a 
considerable displacement of the black hole from the 
center of mass of the binary system: 4.19 ± 0.19 Rq. 
Thus, if the jet axis is aligned with the black hole, 
the base of the jets would orbit about the center of 
mass with a maximum displacement of w 8.4 Rq. A 
highly collimated jet that is being ejected from the 
very near black hole (along the tilted spin axis) as 
the black hole orbits the center of mass would have 
a spiral structure. A more detailed analysis would 
have to be done to see if the wiggles in the radio jets 



with the 3^0 ± 0^2 periodicity observed by Hjellming 
& Rupen 1995 can be entirely explained simply by 
the orbital motion of the black hole. In particular, 
one would have to consider the physical sizes of the 
observed wiggles as well as the complicated physical 
processes taking place inside the jets themselves. 

If the tilt in the jet axis comes about because the 
spin axis of the black hole is tilted, then one needs 
to explain the misalignment of the spin axis of the 
black hole and that of the orbital plane. The Monte 
Carlo simulations of supernovae kicks on neutron star 
binaries by Brandt & Podsiadlowski (1995) predict 
that most of the massive stars in these neutron star 
binaries would have inclined spin axes. If the work 
of Brandt & Podsiadlowski (1995) can be extended 
to the black hole systems, it may turn out that a 
15.5 deg tilt in the black hole's spin axis is not at all 
unusual. 

6.3. Evolution of the System 

Despite the pleasing agreement between single star 
evolutionary models and the observed parameters 
of the secondary, the evolutionary history of GRO 
J1655-40 is clearly much more complex than the evo- 



21 



lutionary history of a single star (e.g. Webbink 1992; 
van den Heuvel 1992; Webbink & Kalogera 1994). 
Webbink & Kalogera (1994) give several conditions 
for the successful formation of a low mass X-ray bi- 
nary, some of the conditions being: (i) the progenitor 
binary must have an extreme mass ratio in order to 
drive the system to a common envelope and the initial 
binary must be wide enough to survive the common 
envelope stage; (ii) the post-common-envelope binary 
must be big enough to accommodate the secondary 
star within its Roche lobe and wide enough to allow 
the primary to evolve to core collapse; and (iii) the 
binary system must survive the core collapse. If the 
conditions stated by Webbink & Kalogera (1994) are 
met and we are left with a 7Mq black hole with a 
2.3 Mq main sequence companion star in an orbit 
with a period of about 2.6 days, the subsequent bi- 
nary evolution will dominated by the evolution of the 
secondary star since the timescale for the evolution of 
that 2.3 Mq star is much shorter than the timescale 
for the loss of orbital angular momentum (King, Kolb, 
& Burderi 1996). The orbit is not likely to be left cir- 
cular after the core collapse, but tidal dissipation will 
circularize the orbit before the secondary can fill its 
Roche lobe (Webbink & Kalogera 1994). Once the 
orbit is circularized, the stellar evolution computa- 
tions we performed show that on a timescale of order 
690 million years, the companion star will evolve off 
the main sequence and attain a radius, luminosity, 
and temperature comparable to the secondary star of 
GRO J 1655-40 that we see today. The radius the star 
attains is roughly the effective radius of the Roche 
lobe of the secondary star in GRO J1655-40. 

GRO J 1655-40 stands out among the black hole 
binaries in that it has a large space velocity (i.e. it 
has a large 7 velocity, see Table [| and Brandt, Pod- 
siadlowski, & Sigurdsson [1995, hereafter BPS95]). 
BPS95 argued that GRO J1655-40 may have ac- 
quired its velocity as the result of a kick caused by 
an asymmetry during the initial collapse of the com- 
pact object (as is thought to be the case with neu- 
tron stars — Brandt & Podsiadlowski 1995; Johnston 
1996). BPS95 proposed that the black hole in GRO 
J1655-40 formed via an intermediate neutron stage, 
and was then converted into a black hole by additional 
accretion from the secondary star or through some 
kind of phase transition in the cooling compact object. 
This scenario predicts that the final black hole binary 
system should have a large mass ratio and that the 
masses of the components would be relatively small 



(Mi ps 3.6 Mq and M 2 ps 0.3 Mq). However, our fits 
to the light curves show that almost the opposite sit- 
uation is true: the mass ratio Q is relatively close to 
unity and the masses of the components are relatively 
large (Mi = 7M and M 2 = 2.3 Mq). Therefore, if 
the compact object passed through an intermediate 
neutron stage, it would have had to accrete ;> 5 Mq 
to attain its present day mass, with the accreted mat- 
ter presumably coming from the companion star. As 
we showed above, the position of the secondary star 
on the Hertzsprung- Russell diagram is consistent with 
the track of a normal 2.3 Mq star halfway between 
the main sequence and the giant branch. If the sec- 
ondary star was once much more massive and it gave 
up most of its mass to the compact primary, its evo- 
lutionary path on the Hertzsprung-Russell diagram 
could be quite different than that of a normal 2.3 Mq 
star. Thus, BPS95's proposed scenario of the for- 
mation of the black hole in GRO J1655-40 seems to 
be less viable in view of the large amount of matter 
the compact object must accrete and the secondary 
star must give up in order to attain their present-day 
masses. If the compact object initially formed as a 
neutron star, it seems likely that a significant amount 
of material from the supernova itself must have fallen 
back, thus converting the neutron star into a black 
hole shortly after its creation. Note that the same su- 
pernova kick invoked by BPS95 to explain the large 
velocity of GRO J1655-40 might also give rise to a 
substantial tilt in the spin axis of the compact object 
as discussed in Section 6.2. 

6.4. Outburst Mechanism 

Recently, King et al. (1996) and van Paradijs 
(1996) discussed the transient behavior in some low 
mass X-ray binaries. In general, a system will be tran- 
sient if the average mass transfer rate M is smaller 
than some critical value. In the case of GRO J1655- 
40 where the main sequence lifetime of the secondary 
star is much shorter than the timescale of the shrink- 
age of the orbit due to angular momentum loss, the 
average mass loss rate is given by King et al. (1996) 
as 

- Mi ps 4 x \{r w P^ m Ml A7 Mq yr" 1 (2) 

where Pd is the orbital period in units of days and 
where the units of M 2 are solar masses. For GRO 
J1655-40, M 2 = 3.4 x 1O" 9 M yr" 1 = 2.16 x 
10 17 g s _1 . This average mass transfer rate is much 
larger than the average mass transfer rates found 
in the other six transient black hole systems, which 
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all have rates around 10 -10 Mq yr -1 (van Paradijs 
1996). For X-ray heated accretion disks, the critical 
mass transfer rate is given by 

M crit « 5 x l(r u M? /3 i£ /8 M© yr" 1 (3) 

where the units of M\ are solar masses and where 
P 3 = P/(3hr) (King et al. 1996). For GRO J1655-40, 
M crit = 1.1 x 1(T 8 Mq yr" 1 . Thus GRO J1655-40 
is not expected to be a persistent X-ray source since 
M < M C rit- It is, however, interesting to note how 
close GRO J1655-40 is to being a persistent X-ray 
source, van Paradijs (1996) gives the following rela- 
tion dividing systems with stable and unstable mass 
transfer: 

log 4, = 35.8 + 1.07 log P(hr) (4) 

where L x is the time-averaged X-ray luminosity over 
one outburst cycle. If a system falls below the line de- 
fined by Equation [| in the P — L x plane, then it will 
be a transient system. If the energy generation rate is 
0.2c 2 per gram of accreted matter (van Paradijs 1996), 
the averaged mass transfer rate of = 2.16 x 10 17 g s _1 
for GRO J1655-40 corresponds to an average X-ray lu- 
minosity of L x = 3.88 x 10 38 erg s _1 . This is slightly 
smaller than the value of L x = 5.31 x 10 38 erg s -1 
predicted from the relation given by Equation (|J). 
The other six transient black hole systems have av- 
erage X-ray luminosities that are at least a factor of 
ten less than the critical X-ray luminosity defined by 
Equation H (see Figure 2 of van Paradijs 1996), which 
suggests that GRO J1655-40 is likely to have more 
frequent X-ray outbursts than the other six transient 
black hole systems. Indeed, there were several ad- 
ditional hard X-ray outbursts after the initial hard 
X-ray outburst observed July, 1994 (Harmon et al. 
1995a), and the soft X-ray outburst observed in late 
April, 1996 (Remillard et al. 1996a) came less than 
one year after the sequence of hard X-ray outbursts 
ended. 

7. Summary 

Using our database of GRO J1655-40 spectra, we 
have established much improved values of the orbital 
period (P = 2^62157 ± 0^00015), the radial velocity 
semiamplitude (K 2 = 228.2 ± 2.2 km s" 1 ), the mass 
function (f(M) = 3.24 ± 0.09 Mq), and MK spectral 
type (F3IV to F6IV). Our photometry taken while the 
system was in true X-ray quiescence shows that light 
from the distorted secondary star dominates, making 



it possible to model the light curves in detail. Our 
model of a distorted secondary star plus a circular ac- 
cretion disk provide excellent fits to the light curves 
taken during true X-ray quiescence. We can add the 
effects of X-ray heating on the secondary star and fit 
the general shape of the March 18-25, 1995 V and / 
light curves, taken during a period of intense activity 
in hard X-rays. The best values of the orbital inclina- 
tion angle i and the mass ratio Q (i = 69.50±0.08 deg 
and Q = 2.99 ± 0.08) may be combined with the mass 
function to give, for the first time, a reliable mass for 
a black hole: M x = 7.02 ± 0.22 Mq. The position 
of the secondary star on the Hertzsprung- Russell dia- 
gram is well determined, and its luminosity, tempera- 
ture, and radius are all consistent with a 2.3 Mq star 
ss 690 million years into its evolution past the ZAMS. 
The average mass accretion rate of the system is much 
larger than the averaged accretion rates of the other 
transient black hole systems, putting GRO J 1655-40 
much closer to the threshold where it would be a per- 
sistently strong X-ray source. 
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A. Appendix: The Eclipsing Light Curve 
Code 

In this Appendix we describe in detail the code 
used to model the eclipsing light curve of GRO J1655- 
40. This program is a modified version of code first 
written by Yoram Avni (Avni & Bahcall 1975; Avni 
1978; see also McClintock & Remillard 1990). We 
will outline here in detail the basic input physics and 
approximations used, not to take credit for the work 
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of Avni and others, but to tell the readers exactly 
what the code does and how it does it. 

A.l. The Potential 

Consider a binary system consisting of a visible 
star of mass Mi and a compact object of mass Mi, 
where the visible star orbits in a circular orbit with 
a Keplerian angular velocity u>f-- Assume the visi- 
ble star is also rotating with an angular velocity loi. 
Following the notation in Avni (1978), we define a 
rectangular coordinate system with the origin at the 
center of the visible star, and with the A-axis point- 
ing towards the compact object, and the Z-axis in 
the direction of Uk- This coordinate system rotates 
with the visible star. The potential can be written as 
(Avni 1978) 



*(x,y,z) = - 



GM 2 



1 m Q n 
1 Qx 

T2 ri 



D 

n 2 (x 2 +y 2 )(l + Q) 



(Al) 



where D is the separation of the two stars, Q = 
Mi /Ma, £1 = LOijijj^, n and r 2 are the distance to 
the centers of the two stars in units of D, and where 
x and y, are X and Y in units of D. If £2 = 1 (the 
visible star is in synchronous rotation), the poten- 
tial given by Equation (Al) reduces to the standard 
Roche potential, and the star can be in hydrostatic 
equilibrium in the rotating frame (Avni 1978). The 
degree to which the star fills its Roche lobe must be 
specified — it is usually taken to be 100%. 

In practice, it is convenient to adopt units of mass 
and distance such that the GM%/D = 1. Typically, 
the star is assumed to be in synchronous rotation, 
so that £2 = 1, and it is assumed that it completely 
fills its Roche lobe. Thus the value of the mass ratio 
Q uniquely determines the function for the potential, 
and hence the geometry of the Roche surface. Once 
the value of Q is specified, the visible star is divided 
into grid points in longitude <f>, where the points 
are spaced equally in cos0, and ANg grid points in 
the co-latitude 9, where the points are spaced equally 
in the angle 6. The value of the potential iff and its 
derivatives are then computed for each point. 

A. 2. Photometric Parameters 

von Zeipel's Theorem (1924) provides a relation- 
ship between the local gravity and the local emergent 
flux in a tidally distorted star. Proofs of von Zeipel's 



Theorem may be found in Kopal & Kitamura (1968) 
and Avni (1978). Using von Zeipel's Theorem, one 
can show that the light emitted from every point on 
the photosphere of the star is the same as the light 
emitted from a plane-parallel atmosphere character- 
ized by the local values of the temperature T e and 
gravity g: 

Tt cc g. (A2) 

As a consequence of von Zeipel's Theorem, the tem- 
perature at any point on the star is given by 



T(x,y,z) 



T, 



pole 



^polc 



p 



(A3) 



where Xpoic and g po ie are the temperature and gravity 
at the pole of the star (i.e. the point on the surface 
of the star where the positive Z-axis emerges) . The 
"gravity darkening exponent" f3 has two values: 0.25 
for stars with radiative atmospheres as shown by von 
Zeipel (1924), and 0.08 for stars with fully convective 
envelopes (Lucy 1967). Thus, to specify the temper- 
ature T(x, y, z) at every point on the star, one must 
input the value of Xp i c . This input temperature is 
usually taken to be the effective temperature of a field 
star with a similar spectral type as the star to be mod- 
elled. In our model, the surface gravity on each point 
on the star can be computed from the derivatives of 
the potential: 



g{x,y,z) = - 



d*ff 
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dy 



dz 



1 1/2 



(A4) 



e.g. (Zhang, Robinson, & Stover 1986). 

Once the temperature is computed for each sur- 
face point, Planck's function is used to approximate 
the relationship between the temperature and the 
monochromatic intensity 



/(A) oc 



exp 



/ he 



\k\T e 



- 1 



(A5) 



where h, c, and k are the usual physical constants. 
The star appears darker near its limb, a phenomenon 
referred to as limb darkening. The code presently uses 
a standard linear limb darkening law expressed as 



/(A) oc 1 — u(X) + u(X) cos ii, 



(A6) 



where fi is the angle between the surface normal and 
the line of sight. The values of the coefficient it(A) 
are taken from standard tables computed from model 
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atmospheres (e.g. Al-Naimiy 1978). The reader is re- 
ferred to Kopal (1959), Kopal & Kitamura (1968), 
and Avni (1978) for more discussions on these ap- 
proximations. 

A. 3. Integration of the Flux from the Star 

Let L(X) be the radiation emitted by the star at 
the wavelength A as seen at a great distance. If 
I(X,x,y,z) is the intensity of the light at a surface 
point, the total observed flux is given by 



(A7) 



1 



= / I(\ x, y, z) cos 7 ds, 



where 7 is the angle of foreshortening, ds is the sur- 
face element, and where the integration is to be done 
over the entire visible surface of the star (Kopal & 
Kitamura 1968). 

To carry out the numerical integration of the flux, 
some quantities need to be defined. First, each point 
on the surface has direction cosines: 



cos cp 
sin <f> cos ( 
sin d> sin 6 



(A8) 



The element of the surface area dS(x, y, z) is given by 

R 2 A<j> Ad 



dS(x,y,z) 



a(x,y,z) 



(A9) 



where R 2 = x 2 + y 2 + z 2 , A<fi and AO are the grid 
sizes in longitude and co-latitude, respectively, and 
where 
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Next, let O be the orbital phase of the observation 
(where = corresponds to the time of the clos- 
est approach of the visible star) and let i be the in- 
clination of the orbit (i — 90deg for an orbit seen 
edge-on). The foreshortening T(x,y,z) of a partic- 
ular point on the star depends on the phase of the 
observation, the inclination, and on the (x, y, z) coor- 
dinates of the point: 
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If the "projection factor" T(x,y,z) < 0, then that 
particular point is not visible. At each phase, the 
flux elements from the visible points (i.e. those with 
T(x,y,z) > 0) are simply summed up: 

L star (e,A)=^]T/(A, x,y,z)T(x,y,z)dS(x,y,z), 
i=i j=i 

(A12) 

where each pair of indices are associated with 
a specific (x,y,z) point on the star. In practice, the 
above sum in Equation ( A12) converges effectively for 
Nfj, > 40 and N e > 14. 

So far, Equations (Al) through ( A12j ) describe the 
basic input physics and mathematics of the original 
Avni code (Avni & Bahcall 1975; Avni 1978). This 
code models the light curve due to a single Roche 
lobe filling star. Extra sources of light such as light 
due to X-ray heating effects are not taken into ac- 
count in the Avni code. To summarize the model so 
far, the user specifies the degree to which the star 
fills its Roche lobe (usually 100%), the value of the 
mass ratio Q, and the rate of the star's rotation (usu- 
ally Q = 1 for synchronous rotation). These three 
quantities define the shape of the surface of the star. 
Then the user must specify the polar temperature of 
the star T po i c , the gravity darkening exponent (3 (ei- 
ther 0.25 or 0.08 depending on whether the star has 
a radiative envelope or a convective envelope), and 
the linearized limb darkening coefficient u(A) appro- 
priate for the star in question in order to compute the 
temperatures over the surface of the star. Once the 
temperatures are known, the Planck function and a 
linearized limb darkening law are used to compute the 
intensities over the surface. Finally, after the orbital 
phase and the inclination of the orbit are specified, 
the total observed flux can be computed. 

A. 4. The Addition of an Accretion Disk 

To make the basic Avni code more realistic, we 
added an accretion disk to the code. Following Zhang 
et al. (1986), the disk is flattened cylinder centered on 
the (invisible) compact object. The plane of the or- 
bit bisects the disk in the ^-direction. The disk is as- 
sumed to be completely optically thick, so that we can 
only see light from its surface and so that anything 
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behind the disk is completely eclipsed. The face of the 
disk is divided up into several grids equally spaced in 
the polar coordinates (r, a). The outer radius of the 
disk rd is scaled to a user specified fraction of the 
effective Roche lobe radius of the primary and the in- 
ner radius r, is set to a very small number (typically 
0.005). Eggleton's (1983) approximation is used to 
compute the effective Roche lobe radius. The thick- 
ness of the disk at the outer edge is Zd = 2rd tan /3 r i m . 
With the exception of user-defined "hot-spots", the 
temperature of the disk does not vary with the az- 
imuth angle. The radial distribution of the tempera- 
ture across the face of the disk is given by 

T(r)=T disk (r/r d )t, (A13) 

where Tdisk is the user specified temperature at the 
outer edge of the disk. For a steady-state, optically 
thick, viscous accretion disk, the power law exponent 
£ = -0.75 (Pringle 1981). As before, Planck's func- 
tion is used to approximate the relationship between 
the temperature of a point on the surface of the disk 
and the monochromatic intensity /disk (A, r, a) of that 
point. 

Recently Diaz, Wade, & Hubeny (1996) discussed 
the importance of including corrections for limb dark- 
ening in models of disk spectra. Limb darkening is im- 
portant for disks since the limb darkening corrections 
depend on the temperature (and hence the radius in 
our model) and on the effective wavelength. There- 
fore the slope of a blackbody disk's spectrum will be 
in error if no limb darkening corrections are made. 
An important result of the work of Diaz et al. (1996) 
is that the limb darkening law in the optically thick 
rings of their accretion disk models is very similar to 
the limb darkening law in a stellar atmosphere. As a 
result, one may use the same linearized limb darken- 
ing coefficients for disks that one uses for stars. We 
therefore have 

7disk(A,r,a) cx 1 - u dis k(A, T(r)) 

+u disk (X,T(r))cosi (A14) 

where the values of Udisk(A, T(r)) are interpolated 
from tables given by Wade & Rucinski (1985). 

The foreshortening angle of the normal of a surface 
element on the face of the disk is simply cosi. In the 
absence of eclipses, the observed flux from the face of 



the accretion disk is: 

L d isk(B, A) = ^^idi S k(A,r,a)(cosi)rAr Aa, 

i=i j=i 

(A15) 

where N r and N a are the number of grid points in ra- 
dius and azimuth, respectively, where Ar and Aa are 
the grid spacings in radius and azimuth, and where 
each pair of (i, j) indices are associated with a specific 
(r, a) point on the disk. 

To handle the situations where the disk has a sub- 
stantial thickness, light from the rim is also accounted 
for. The rim is divided into 11 grid points in the z 
direction, and N a grid points in azimuth. The tem- 
perature of the rim is Td; s k, and the intensities of each 
rim point are found using Planck's function. If the az- 
imuth angle a is measured from the A-axis, then the 
foreshortening factor r r i m (a) of a point on the rim at 
the orbital phase is 

rrim(a) = sini cos(a + 9). (A16) 

If r r i m (o;) < 0, the point is hidden. A limb darkening 
correction is made for the edge surface elements based 
on the value of r r j m (a): 

irim(A, a) CX 1 - Mdisk(A,T d i s k) 

+«disk(A,Tdi s k)r rim 

(a) (A17) 

The total flux from the rim is found by summing the 
intensities of all visible points: 

11 N a 

rim(tt) Az Aa, 

i=l j=l 

(A18) 

where Az is the step size in the Z direction, and 
where each pair of indices are associated with a 
specific (z, a) point on the rim. 

Hot spots can be added to the disk. The spots are 
confined to radii r cut < r < r d and azimuth angles 
ot\ < ex. < a\, where at the moment i — (no spots), 
1 or 2. The temperature of the spot T spot is added to 
the temperature T(r, a) of each point on the disk rim 
and face that lies inside the spot boundaries. The in- 
tensities of each point are computed as before and the 
i nteg rations are carried out as before (e.g. Equations 
( [All ) and (|AL§)). 

In the absence of eclipses, the total observed flux 
from the star and the disk is 

itotai(©, A) = L star (e, A) + i dis k(e, A) + L rim (6, A). 

(A19) 
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A. 5. The Addition of Eclipses 

In the case where the inclination angle i is large 
enough to allow for eclipses, the procedure for com- 
puting the light curves must be modified slightly. 
First, the intensities of each point on the star, disk 
face, and disk rim are computed from the list of in- 
put parameters. We can determine which body is in 
front of the other body from the orbital phase — 
the star will obviously be closer than the disk for 
Odeg < 6 < 90deg or for 270 deg < 6 < 360 deg 
and the other way around for 90 deg < 8 < 270 deg. 
The integration of the flux elements over the visible 
surface of the object in front is carried out as before. 
Then, all points on the edge of the body in front (i.e. 
its horizon) are determined and are projected onto the 
plane of the sky. The horizon of the eclipsing body 
defines some polygon on the plane of the sky. Then, 
as the integration of the flux elements of the body in 
back is carried out, each potentially visible point is 
projected onto the plane of the sky and checked to 
see if it falls inside the horizon of the eclipsing body. 
If the point is eclipsed, it is not included in the inte- 
gration. Once all of the visible points are included in 
the integrations, t he to tal observed flux is found as 
before (Equation ( A19 )). 



A. 6. The Addition of X-ray Heating 

If the secondary star intercepts a large amount of 
flux from the other body (either X-ray heating from 
the vicinity of the compact object or optical light 
from the disk), there may be a large amount of re- 
processed radiation. The irradiating flux can strongly 
alter the distribution of temperatures across the face 
of the star. A rigorous treatment of these "reflection 
effects" requires difficult and time-consuming compu- 
tations (e.g. Kopal 1959; Wilson & Devinney 1971). 
Following Zhang ct al. (1986), we will use a simpli- 
fied treatment of the heating effects. We will consider 
only heating of the star by X-ray heating, and not the 
heating of the disk by the star. For simplicity, all of 
the irradiating flux Fj„ due to the X-ray heating is 
assumed to come from the center of the compact ob- 
ject. 

The total energy in X-rays intercepted by the 
Roche lobe filling secondary star does not depend on 
the scale of the orbit since the surface area of Roche 
lobe filling star facing the X-ray source is proportional 
to the square of the orbital separation. However, the 
change in the star's atmosphere caused by the X-ray 



heating does depends on the size of the orbit. As the 
distance between the star and the X-ray source is in- 
creased, each square centimeter on the star's surface 
receives less energy in X-rays. The change in the lo- 
cal temperature depends on the amount of absorbed 
X-ray flux per unit area. So to find the magnitude 
of the irradiating flux seen by each point on the sec- 
ondary, the X-ray luminosity of the X-ray source L x 
(in ergs s _1 ) and the physical size of the orbit must 
be specified. The code uses the mass function f(M) 
and the orbital period P to find the scale of the orbit. 

The distance d (in units of the orbital separation) 
between each point on the face of the star and the 
X-ray source can be found 



d(x, y, z) = (1 + z 2 + R 2 - 2R£ X 



(A20) 



where i x is one of the direction cosines (Equation 



(AS)), and where R = x + y + z 2 . The foreshort- 



ening factor Tx-ra,y(x, y, z) is given by 

x(l — x) — y 2 + z 2 



r X -T ay (x,y,z) 



Rd 



(A21) 



where Tx-ia,y(x,y, z) < indicates that the point on 
the star cannot see the X-ray source. The amount of 
X-ray flux each surface element intercepts is 



F irl (x,y, z) = 



L x 



Tx-Ta.y(x,y,z). 



4ira 2 d(x,y,z) 2 J * A ~ rayV 

(A22) 

The modified temperature of each surface point that 
can see the X-ray source is given by 



^X-rayC 2 -' 2/' Z ) 



pole 



9 P o\c{x,y,z) 



if) 



WF iTI 



a 

(A23) 

where W is the albedo, and a is the Stefan-Boltzmann 
constant (Zhang et al. 1986). The albedo W is not 
very well determined, and it is usually set to 0.5 (see 
Zhang et al. 1986). 

In the cases there the disk has a substantial thick- 
ness, some points on the face of the star may be be- 
low the rim of the disk as seen from the compact ob- 
ject. These points would not receive any X-ray flux 
from the compact object and hence their tempera- 
tures would not change. It is easy to find the angle 
above the plane that a point on the face of the star 
has (as seen from the compact object): 



&(x, y, z) — arctan 



(1 — x) cosrj 



(A24) 
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where r\ = arctan(y/x). This angle <j>(x,y, z) is com- 
pared to the angle of the rim /3 r i m to see whether 
or not the point on the star is blocked by the disk. 
If /3 r im > &(x,y, z), the point cannot see the X-ray 
source. 

Once each point on the face of the star has been 
assigned a new temperature based on the amount of 
X-ray flux it sees, the intensities and the fluxes are 
found following the procedures outlined above. The 
net result is that the points on the star facing the X- 
ray source are brighter than they would otherwise be 
because the temperatures were raised. 

A. 7. Filter Response Curves 

We compute the total flux L tota \(Q, X) at each 
phase for several values of the effective wavelength 
A using limb darkening coefficients w(A) interpolated 
from tables given by Wade & Rucinski (1985). This 
model "spectrum" is integrated with the standard 
UBVRI filter response curves (Bessell 1990) to pro- 
duce model UBVRI magnitudes. In this way we 
can perform simultaneous fits to light curves in sev- 
eral colors. In the case of GRO J1655-40, there is a 
large amount of reddening in the light curves. Rather 
than attempting to include the reddening correction 
in the fits, we instead fit the models to the normalized 
BVRI data. In this way, the error of E(B — V) does 
not enter into the analysis. We do, however, lose the 
absolute color information contained in the model. 

A. 8. Light Curve Fitting and Error Analysis 

We first compute a grid of models and x 2 fits to the 
data using a wide range of parameters spaced in uni- 
form steps. Then several sets of parameters that pro- 
duce low x 2 values are chosen and are used as start- 
ing points for an optimization routine based on the 
GRIDLS program given in Bevington (1969). This 
"grid search" method was picked because it is easy to 
implement and it does not require the computation 
of derivatives. After several runs with different start- 
ing points, the set of parameters that gives the lowest 
overall \ 2 fit to the data is chosen. By using several 
different starting points for the optimization routine, 
we can be reasonably sure that we have found the 
global minimum of \ 2 , rather than a local minimum. 

We use three methods to estimate the statistical er- 
rors on the fitted parameters. The first one is a "boot- 
strap" method (Press et al. 1992). Artificial data sets 
are constructed from the original data set and the ex- 



act fitting procedures which were used on the original 
data are used on the artificial data. The la errors 
on the parameters are then derived from the distribu- 
tions of the fitted parameters from the many artificial 
data sets. The second technique we use is outlined by 
Zhang et al. (1986). Basically the la error on a fit- 
ted parameter ao is found by fixing the value of that 
parameter at a new value a' = a + S and finding 
the values of the rest of the parameters a\ through 
a n that minimize \ 2 - The probable error of ao is re- 
lated to the change in \ 2 an d to 5. We typically use 
two or three different values of 5 (1%, 2%, and 3% of 
the value of the parameter) and average the resulting 
values of a. Finally, the GRIDLS program given in 
Bevington (1969) also provides error estimates based 
on the shape of the chi-square hypersurface near the 
minimum. We used primarily the errors derived from 
the bootstrap method. We stress here that the er- 
rors estimated are only internal statistical errors, and 
not systematic errors which we have no easy way to 
compute. 
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